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ABSTRACT 

We report high-resolution optical spectroscopy of the low-mass X-ray binary V395 Car/2S 0921-630 ob- 
tained with the MIKE echelle spectrograph on the Magellan-Clay telescope. Our spectra are obtained near 
inferior conjunction of the mass donor star and we exploit the absorption lines originating from the back- 
side of the K-type object to accurately derive its rotational velocity. Using KO-Kl 111 templates, we find 
vsin/ = 32.9 ±0.8 km s"'. We show that the choice of template star and the assumed limb darkening coefficient 
has little impact on the derived rotational velocity. This value is a significant revision downwards compared 
to previously published values. We derive new system parameter constraints in the light of our much lower 
rotational velocity. We find Mi = 1 .44 ± 0. lOM©, Mi = 0.35 ± 0.03Mq, and q = 0.24 ± 0.02 where the errors 
have been estimated through a Monte-Carlo simulation. A possible remaining systematic effect is the fact 
that we may be over-estimating the orbital velocity of the mass donor due to irradiation effects. However, 
any correction for this effect will only reduce the compact object mass further, down to a minimum mass of 
Ml = 1.05 ±0.08Mq. There is thus strong evidence that the compact object in this binary is a neutron star of 
rather typical mass and that the previously reported mass values of 2-4M0 were too high due to an over-estimate 
of the rotational broadening. 
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1. INTRODUCTION 

Low-mass X-ray binaries (LMXBs) are binary systems in 
which a neutron star or a black hole accretes matter from a 
low-mass companion star These systems provide a labora- 
tory to test the behavior of matter under physical conditions 
that are unattainable on Earth. One of the ultimate goals of 
the study of neutron stars is to determine the equation of state 
(EoS) that describes the relation between pressure and density 
of ma tter under the extrem e conditions encountered in neutron 
stars dLattimer & Prakashi 2001 , 2004). As a result of the cur- 
rent paucity of observational constraints there are many the- 
ories of the EoS of matter at neutron-star densities. Through 
the measurement of the masses and radii of neutron stars these 
theories can be tested. For a specific EoS, one can derive a 
firm upper limit on the mass of the neutron star, above which 
the object is not stable and would collapse into a black hole. 
Neutron stars with masses well above 1 AMq cannot exist for 
so-called soft EoSs. Therefore, measuring a high mass for 
even one neutron star would impl y the firm rejection of many 
propo sed EoSs (see discussion bv lvan Paradiis & McCUntocH 
Il995h . 

V395 Car/2S 0921-630 is a promising accreting binary sys- 
tem for accurate system parameter work . It shows regular 
eclipses in the X-rays (Mason et a l.lll987|) as well a s dips in 
the optical lightcurve (Chevalier & Ilovaiskylll98Th . These 
reveal a 9.02d orbital period and imply a high orbital incli- 
nation. The unknown orbital inclination is often a signifi- 
cant handicap when trying to determine binary parameters, 
but in V395 Car the sin; factor is well constrained. In the 
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opti cal, the K-type donor star absorpti on fine spectrum is vis- 
ible (iBranduardi-Raymont et al.lll983h which allows us to de - 
termine its rotational velocity, vsin/ ('Shahbazet al. 1999lX 
as well as the radial ve locity amplitude K2 (Shahbaz et al] 
2004; Jonker et al.ll2005l) . These previous studies suggested 
the presence of a rather massive compact object with a mass 
between ~2-4M0. In the case of V395 Car the main source 
of error in the mass of the compact object is the large error on 
the rotational velocity (Jonker et al. 2005) which was de rived 
from relatively low resolution data (iShahbaz et al.lll999h . 

Here, we present Magellan Inamori Kyocera Echelle 
(MIKE) data of the V^16 magnitude counterpart of V395 
Car/2S 0921-630 to determine the rotational velocity of the 
companion star accurately. We show in this Letter that we 
find a significantly lower value for the rotational broadening 
than previously determined, and explore the effects on the in- 
ferred neutron star mass. 

2. OBSERVATIONS & REDUCTION 

We observed V395 Car/2S 0921-630 with the MIKE 
echelle spectrograph mounted on the Magellan Clay tele- 
scope at Las Campanas Observatory. Four exposures of 30 
mins each were obtained on Jan. 25, 2005 (MJD 53395) 
between 7:05-9:15 UTC. We employed the spectrograph in 
its dual-beam mode using a 1" wide slit and a dichroic. 
The 2048x4096 pixel CCD detectors were binned on-chip in 
2x2 mode. This delivers a wavelength coverage of 3430- 
5140A on the blue MIT detector with a spectral dispersion 
of 0.04A/pixel, while the red SITe detector covered 5220- 
9400A at 0.1 OA/pixel. The spectrophotometric standard HR 
4468 was observed immediately after the target exposures. 
Observing conditions were good with clear skies and 0.5" see- 
ing which resulted in a seeing limited resolution of ^0.08 A in 
the blue and 0.17 A in the red. Exposures of 5 bright template 
stars of spectral type K were obtained with the same setup a 
few nights earlier on Jan. 21 2005. 
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TABLE 1 

Rotational broadening as a function of template 



Star 


Spectral 


I' sin ! (km s ' )" 




ID 


type 


0.5 


0.75 


HD 100623 


KOV 


30.5±0.8 


31.5±0.9 


HD 124106 


KIV 


30.5±0.8 


31.4±0.8 


HD 99322 


KOIII 


32.3±0.7 


33.6±0.8 


HD 121416 


Kim 


31.3±0.7 


32.5±0.8 


HD 83240 


KlIII 


31.4±0.8 


32.7±0.9 


MS mean 




30.5±0.8 


31.5±0.9 


Giant mean 




31.7±0.7 


32.9±0.8 



The data were reduced using the MIKE REDUX^ software 
package that is tailored for the instrument. It performs the 
standard CCD calibration steps, constructs a wavelength so- 
lution using ThAr exposures obtained at the position of the 
target and optimally extracts the object frames. The orders 
were merged using our spectrophotometric flux standard and 
then exported into our analysis software (MOLLY). We re- 
binned all our object and template star spectra onto the same 
wavelength scale using a uniform velocity scale and their con- 
tinua were normalised using spline fits while masking strong 
spectral features. All spectra were shifted into a heliocentric 
velocity frame by removing the Earth's velocity. 

3. ANALYSIS 

The principle aim of our observations was to accurately de- 
termine the rotational velocity of the mass donor star. We 
intentionally obtained the observations close to inferior con- 
junction of the mass donor star when its contribution to the 
overall light is expected to be the largest. Another advantage 
is that at this phase we observe the back side of the donor 
star. This side is not subjected to X-ray heating due to the 
luminous accretion flow around the compact object and the 
non-spherical effects due to its Roche-lobe shape are mini- 
mal. According to the ephemeris of Jonker et al. (2005), our 
observations span 0.89-0.90 in orbital phase. 

Any intrinsic velocity shifts between our target spectra and 
the template stars were first removed by cross-correlating the 
normalised spectra while masking the emission lines in our 
target spectra as well as regions containing telluric absorption 
features. The template stars were then rotationally broadened 
by a varying amount using a spherical model and including 
a linear limb darkening coefficient. We then subtracted these 
broadened templates from our individual target spectra while 
optimising for a variable scaling factor / reflecting the fact 
that not all the light in the system originates from the mass 
donor. This optimal subtraction is performed by minimizing 
the root mean square (rms) residuals of the data - (/x constant 
+ (l-/)x template spectrum). Again the regions containing 
telluric absorption and H and He emission lines in the object 
spectrum were masked. The goodness of fit of this subtraction 
is then evaluated as a function of the applied rotational broad- 
ening using the x^-statistic. We plot two small segments of 
our observed average spectrum of V395 Car together with a 
broadened template spectrum in Figure [T] 

We paired up the individual spectra with our 5 templates 
and found that in all cases, the optimal rotational broadening 
was found to be close to 32 km s"', a value that is signifi- 
cantly lower than the 64 ± 9 km s"' reported by Shahbaz et al. 
(1999). We found no significant differences between the val- 
ues derived from the individual exposures and thus averaged 

^ http://web.mit.edu/~burles/www/MlKE/ 
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Fig. 1 . — Two small segments of data showing the mean spectrum of V395 
Car together with a broadened KOIII template. A large number of resolved 
absorption lines that are well matched by early K-stars are clearly visible in 
both wavelength regions. The bottom panel is from the blue arm data and also 
contains the strong Hell 4686 emission Une as well as Hel 4713 absorption 
originating in the accretion flow . 

our four shifted spectra together using variance weights to 
boost the S/N. We also found that the red spectra are preferred 
since they contain wide spectral windows free from emission 
lines where strong K-star absorption features are expected. 
For our final analysis, we therefore repeated our rotational 
broadening analysis using the mean red spectrum of V395 Car 
and broadening all templates between 25-50 km s"' in steps 
of 0.5 km s"'. The optimal vsin/ value for each template and 
for a given choice of the limb darkening coefficient was de- 
rived by fitting a cubic function to the x versus vsinj values 
(Figure 12^). Our best fit solution usually delivered a xi ^ 2 
reflecting the rather simplistic assumption that the spectrum 
can be modeled by a constant and featureless accretion disk 
spectrum and the mass donor star as represented by our spec- 
tral templates. In order to determine more reliable error esti- 
mates on our vsin / values, we used a Monte-Carlo approach. 
We simulated 500 copies of our input target spectrum using a 
bootstrap technique where the input spectrum is resampled by 
randomly selecting data points from it. Each bootstrap copy 
always has the same number of data points as the input spec- 
trum by allowing points to be selected multiple times or not 
at all. For each bootstrap copy, the optimal vsin / is evaluated 
as described above. We found that the vsin / values were well 
described by a Gaussian distribution (Figure|2j)), and thus cal- 
culated the mean and rms of our 500 vsin / values to deliver a 
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Fig. 2. — a) An example showing the dependence of the achieved x?y 
statistic for an optimally subtracted template star as a function of the applied 
rotational broadening vsini. Symbols show the actual values returned from 
our fits, while the solid line is a cubic fit that is used to determine the best 
vsini for each run. b) Histogram of the derived rotational broadening val- 
ues calculated from 500 trial runs of subtracting the KlIII star HD 121416. 
The values are well characterised by a Gaussian distribution as shown by the 
model Gaussian and provides the final v sin i and its error for a given template, 
31.3±0.7 km s"' in this case. 

V sin / value and its 1 a error for each template star. We provide 
the results for our 5 template stars in Table [T] for two values 
of the limb darkening coefficient. For K-giants, one expects 
a typical limb dark ening coefficient of 0.75 in the continuum 
dClaret et al.lll995h . but we also list the equivalent values for 
an assumed limb darkening coefficient of 0.50 in order to de- 
termine its effect on the derived vsin/ values. 

4. IMPLICATIONS 

Thanks to the high spectral resolution and good S/N of our 
Magellan data, we were able to extract accurate v sin / values 
with < 1 km s"' errors. The value was also stable against 
choice of template, although the required vsin / for the giant 
templates is on average slightly higher than that for the main- 
sequence templates (TablefTli. This likely reflects the fact that 
the absorption lines found in the unbroadened main sequence 
templates have higher pressure atmospheres and thus are in- 
trinsically broader and require a slightly smaller amount of 
additional rotational broadening. This difference is < 2(t as 
is the difference between individual templates. It can also be 
seen that the effect of limb darkening is relatively modest, 
shifting the vsin/ values systematically down by 1 km s"' if 
we lower the coefficient down to 0.5. This trend continues if 
we push the line limb darkening further down and in the ex- 
treme lowers vsin / by 3 km s~' (see also Shahbaz & Watson 



2007). Since we expect the Roche-lobe filling donor star in 
its 9 day orbit around the neutron star to be closer to a giant 
rather than a main sequence dwarf, we use vsin/=32.9±0.8 
km s~' in the remainder of this Letter This is the mean value 
derived for our three giant templates while using a limb dark- 
ening of 0.75, as is appropriate for the R-band continuum in 
early K-giants. We remark that the bias introduced by this 
choice is very modest since all individual values are within 
2a of this and we will see below that the uncertainty on the 
derived neutron star mass is no longer dominated by the un- 
certainty in vsin/, but instead by the radial velocity amplitude 
K2. Very similar results were recently obtained by Shahbaz & 
Watson (2007) using VLT echelle spectroscopy. 

As mentioned above, these recent vsin / values are a factor 
of 2 lower than the previously reported value of Shahbaz et 
al. (1999). Shahbaz & Watson (2007) conclude that the dom- 
inant absorption blend near 6495 A may have biased the 1999 
vsin / determination which used data with relatively low sig- 
nal to noise as well as low spectral resolution. In our case, the 
high spectral resolution ensured that neither the instrumental 
profile of the spectrograph nor variable slit illumination ef- 
fects had an impact on the analysis since the rotational broad- 
ening exceeded our resolution by a large amount. We also 
used templates that were obtained with the same instrument 
and had sufficient signal to noise to detect and resolve a large 
number of absorption lines. This permitted the determination 
of reliable values for vsin/ with a better than 1 km s"' statis- 
tical precision. 

4.1. The masses of the stellar components 

Armed with an improved determination of vsin/, we can 
now derive new constraints on the mass of the compact object 
and its K-type donor The measured radial velocity semi- 
amplitude K2 and the rotational velocity give the mass ra- 
tio, q = M2/M1 of th e system via ^ = QA6[(\+qfq\^ 
(I Wade & Ho rne 1988'). We can then solve for the masses of 
the stellar components using the mass functions which deliver 
Msin^ /. Deriving actual masses requires knowledge of the bi- 
nary inclination, /, which for V395 Car is thought to be close 
to 80 degrees (Mason et al. 1987). Values for K2 have been 
reported in Jonker et al. (2005) which found 7^2=99. 1±3.1 
km s-i while Shahbaz et al. (2004) quote /:2=92.9±3.8 km 
s"' . This observed amplitude may be subject to a so-called K- 
correction due to the fact that the inner face of the companion 
star can be heated by X-ray emission coming from near the 
compact object. The center of light measured via the spectral 
absorption features is then offset from the true center of mass 
causing an overestimate of A'2. The ma gnitude of this effect 
depen ds on the irradiation geometry (e.g. iMunoz-Darias et alj 
l2005h as well as the intensity of the incoming radiation. Shah- 
baz et al. (2004) argue that this correction is not significant in 
the case of their K2, but Jonker et al. (2005) provide some 
evidence that X-ray heating does seem to be play a role in 
V395 Car. In the most extreme case of no absorption line 
contribution from the heated front face of the Roche lobe, the 
K-correction can be as large as /S.K < 4/3tt x vsin / = 14 km 
s-i. 

In order to determine the uncertainties on the component 
masses given the observables K2, vsin / and /, we ran a Monte- 
Carlo simulation. The X-ray eclipse discussed in Mason et 
al. (1987) gives us a good handle on the inclination. At face 
value, /r2=99.1 km s"' and vsin/=32.9 km s"' gives a mass ra- 
tio of ^ = 0.25. Using Table 2 in Mason et al. (1987) suggests 
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Fig. 3. — The four panels plot the distribution of system parameters for 
5000 randomly selected combinations of K2 and vsinr. Clockwise starting 
at the upper left corner, we show the mass of the companion star, the mass 
ratio q, the radial velocity amplitude Ki and v sin ; plotted against the neutron 
star mass. The thick dashed curves show in which direction these probability 
clouds move when an irradiation correction to K2 is applied with the solid 
circles indicating the system parameter solutions at the two extremes of no 
correction and maximum K-correction. 

that for this mass ratio, fits to the X-ray eclipse lightcurve 
imply an inclination angle of ~83 degrees. Although some- 
what model dependent, the inclination simply enters as sin^ /, 
which is only a correction at the 2% level for these high in- 
clinations. Using / = 83°, we select random values for K2 
and vsin/ by picking values from a normal distribution cor- 
responding to the mean and la error of the observed values. 
The system parameters are then calculated for each random 
set of parameters, and the process is repeated 5000 times. 
Figure [3] shows the corresponding probability clouds for var- 
ious parameters. In order to illustrate the effect of the un- 
known K2 correction due to X-ray heating, we plot the clouds 



for zero K-correction (/r2=99.1±3.1 km s~') and then show 
with dashed curves the steady trajectory the probability clouds 
would take as a function of the K-correction up to the maxi- 
mally allowed correction. We find that at the canonical values 
for/r2 andysin/,the neutron star mass isMi = 1.44±O.1OM0, 
M2 = 0.35 ± O.O3M0 and q = 0.24 ± 0.02. The uncertainties 
reflect the 1(t/68% error values obtained by propagating the 
la error values on the input parameters through our Monte- 
Carlo simulation. Our revised vsin/ value has thus brought 
the neutron star mass down to a value remarkably close to 
the canonical neutron star mass. From the curves in Figure 
[3] we see that the neutron star mass would go down fur- 
ther if any K-correction is applied with a minimum mass of 
Ml = 1.05 ±0.08Mq at the maximal K-correction. The mass 
for the donor star is also much lower than previously reported. 
The evolutionary timescale for such a low mass star is rather 
long and it is clear that single star evolution is not able to 
produce a 0.35Mq object that is able to fill its Roche-lobe in 
a 9 day orbit. This is not an uncommon finding in LMXBs 
and might be due to a phase of unstable and non-conservative 
mass transfer that caused a significant amount of mass to be 
lost from the donor. 

We conclude that we have derived a revised estimate for 
the rotational broadening of the K-type mass donor star in the 
neutron star binary V395 Car. From this we derive new con- 
straints on the mass of the neutron star which is most likely 
very close to the canonical neutron star mass of 1 AMq. V395 
Car is thus yet another case where a potentially massive neu- 
tron star appears to be brought back into line with the bulk of 
the known neutron stars. 
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